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Abstract 

Radiation damping of the motion of charged particles in relativistic, optically thin plasmas is 
described within the framework of the covariant gyrokinetic theory Q|. It involves description of 
the collisionless single-particle dynamics as well as the Vlasov and Maxwell equations both written 
in the covariant formulation. The damping causes corrections to the phase-space trajectory of 
the particle, as well as to the form of the kinetic equation itself, due to the failure of conditions 
of the Liouville theorem. Both effects result independent of the gyrophase, which is retained as 
an ignorable variable. In addition, the applicability range of the covariant gyrokinetic theory is 
extended to describe short-wavelength perturbations with the background of zero parallel electric 
field. The presented theory is suitable for description of magnetized, relativistic, collisionless 
plasmas in the context of astrophysical or laboratory problems. Non-uniquenes of the gyrokinetic 
representation and consequences thereof are discussed. 
PACS: 52.30.Gz, 52.27.Ny, 04.20. Cv 
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I. INTRODUCTION 

Gyrokinetic description of hot magnetized plasmas P] is an effective tool that allows to 
combine relative simplicity of the drift approximation with adequate account for the finite- 
Larmor-radius (FLR) effects, including description of short-wavelength perturbations. Cur- 

nnnn 

rently, most applications are in the tokamak physics, p, |9|, [lOl, ll8| where it is used for 
of transport, stability and turbulence. Earlier relativistic modifications of the theory[ 
were designed for simulation of confinement of fusion products. Meanwhile, the covariant, 
: our- dimensional formalism of the gyrokinetic (GK) theory, developed in our previous papers 
l|, aims at being capable of describing fast flows of matter, such as taking place in 

accretion disks and relativistic jets in the realm of astrophysics. In particular, it is consistent 
with the general relativity, permits the use of arbitrary four-dimensional coordinate systems 
and non-canonical phase-space variables. Since four- dimensional formulation utilizes inter- 
nal symmetries of electrodynamics, so that the applicability range of the gyrokinetic theory 
extends even in its no n- relativistic limit. 

Four-dimensional representation of drift trajectories, and its relationship to symmetries of 
the Faraday tensor were first published by Fradkin. [16] On that basis, using the Lie-transform 
methods, Boghosian|17] constructed the first version of the covariant, four-dimensional gy- 
rokinetic transformation. Still, that version was not covariant in the general-relativistic 
sense, inherited several restrictions of the non-relativistic theories, and was too complicated 
to be of practical use. Another formulation of covariant gyrokinetic theory was launched 
by Beklemishev and Tessarotto^^ using a perturbative Lagrangian approach introduced 
by Littlejohn.P] The theory was just first-order in the Larmor radius, but it was simple, 
straightforward and fully covariant. An important point of that paper was concerned the 
non-uniqueness of the gyrokinetic transformati on, on which we are going to elaborate fur- 
ther below. Additional work in this direction Q,^^ extended analysis to second order in the 
Larmor radius and in the amplitude of the wave-field, added explicit solutions for particle 
trajectories, determined the covariant gyrokinetic equation, and the effective source-term 
for the Maxwell equations. ^ 

The covariant GK theory, published in has two important limitations, which we aim 
at amending in the present paper. The first is typical to all existing gyrokinetic theories, 
namely, their inability to account for cyclotron losses. The losses, while relatively small in 
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the non-relativistic limit, become significant for relativistic particles. In both cases they 
destroy the most important integral of motion exploited by GK theories, the "magnetic mo- 
ment". The effect, however, is difficult to account for within the context of any Lagrangian 
formulation, since it destroys the Hamiltonian nature of the system. Hence its inclusion in 
the context of GK theory requires a specific treatment. 

Radiation damping of the motion of a charged particle is a natural process described in 



textbooks. 



2l| However, there is a question whether it can be directly applied for plasma 



particles as is. Indeed, the effect is resonant and definitely falls out of the scope of the theory 
for cold dense plasmas, where the cyclotron emission is trapped and becomes a collective 
process. This statement can be illustrated by the following example: a single particle moving 
in a circle emits waves, while a steady circular current (or lots of particles moving in the same 
circle) emits nothing. Well, the opposite limit is also realistic. In hot plasmas relativistic 
effects destroy the resonance, so that the radiation losses grow and lose coherence, i.e., occur 
as if each particle is independent. This is the case we are going to explore: a hot rarefied 
plasma, with a small optical depth in the cyclotron/synchrotron frequency range. The "small 
optical depth" clause is included so that the radiation can escape, and the inverse processes 
of scattering/absorption can be neglected. 

In a hot rarefied plasma, the radiation damping causes corrections to the phase-space 
trajectories of particles, as well as to the form of the kinetic equation itself, since the phase- 
space volume is no longer conserved. We describe these effects as small corrections, evalu- 
ating them in the unperturbed gyrokinetic variables. Both effects result independent of the 
gyrophase, which thus can be retained as an ignorable variable. 

The second limitation of the covariant GK theory developed in , which we 

intend to overcome in this paper, is its inability to treat the short-wavelength perturbations 
when the background parallel electric field is zero. This effectively limits description of typ- 
ical plasmas in fusion devices, where the parallel electric field is small, to long-wavelength 
perturbations. In contrast, all alternative formulations of the gyrokinetic theory are appli- 
cable only if the electric field is small. By searching for the gyrokinetic transformation in 
a broader class than before, we are able to derive an alternative formulation, which is valid 
in all limits, i.e., both for small and for large (of the order of the magnetic field) parallel 
electric fields. 

We also try to look at the problem of alternative formulations of the gyrokinetic theory 
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in broader context: is the gyrokinetic description unique? Are all alternative formulations 
equivalent in some sense in their common applicability range? We intend to show that non- 
equivalent descriptions arc possible and still valid in certain orders of the perturbation theory. 
This feature is not related to the Lagrangian approach but is common to all perturbative 
methods. 

Direct comparison of two non-equivalent gyrokinetic theories, in the common applicability 
range, shows that even such basic feature as drift trajectories of the guiding center are 
affected by the choice of free parameters. Exploration of these degrees of freedom provides 
new insight into the nature and applicability of gyrokinetic theories. For example, a closed 
drift trajectory in one representation will not necessarily be closed in the other one. The 
same applies for drift surfaces in magnetic traps as well. This feature is very important for 
theoretical interpretation of cross-field transport processes. Indeed, it follows that though 
existence of closed drift surfaces means that particles are contained, the reverse is not true 
- those closed surfaces might still exist in other representations of the drift variables. 

The paper is organized as follows. In Section 11 we start by introducing necessary nota- 
tions and definitions, and in Section 111 present the derivation of the covariant gyrokinetic 
transformation, valid for all values of the electric field, and the relevant equations of motion. 
On this basis, in Section IV we discuss the uniqueness and convergence or the gyrokinetic 
approximation. Next, in Section V, we describe the effect of the radiation damping on the 
motion of a particle in the gyrokinetic variables. In Section VI the covariant kinetic and 
gyrokinetic equations arc derived with due account for the radiation damping effects. The 
section also contains our proof of the Liouville theorem in Lagrangian variables, and explicit 
expressions for the rate of contraction of the phase-space volume. Results are summarized 
in Conclusion. 

Throughout the paper the following notations are used: for four-dimensional vectors we 
use bold letters, for tensor fields and operators use roman and calligraphic fonts. Normal 

math letters are either scalars, or components of vector and tensor fields. 

II. COORDINATES, THE TETRAD AND THE VARIATIONAL PRINCIPLE 

In order to carry out the program indicated in the introduction we first introduce the 
phase-space coordinates, and the tetrad of unit four-vectors linked to symmetry planes of 
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the Faraday tensor. In the covariant formulation the tetrad plays a role similar to magnetic 
coordinates. In these terms we formulate the variational principle for description of the 
single particle dynamics. Finally, we introduce and discuss the template and the ordering, 
which will be used later for construction of the gyrokinetic transformation of the phase-space 
variables. 

Covariance of the representation is defined as follows: suppose we have arbitrary, user- 
defined coordinate system, of which we know the metric tensor, g^^,, and the time-like 
invariant interval ds is defined as 

ds^ = g^.dx'^dx'', (1) 

where the Greek indices are assumed to go through /x, = 0,...,3. The signature of g^i, 
is assumed to be < H-,3— >, so that the zero coordinate is also time-like, and s is the 
same as the "proper time" .The four- velocity u is then defined in terms of its countervariant 
components according to 

= HI- 

As a result of (pQ), u is defined as the unit 4- vector tangent to the trajectory. It is time-like, 

= 1. (3) 

Thus, for an arbitrary coordinate system we have defined the seven- dimensional phase- 
space as {x^^Uy). If, then, we present an expression for a particle trajectory in this phase- 
space, it will be covariant, i.e., valid in any coordinates. Section VI contains more informa- 
tion about the phase space and distribution functions in it. 



A. The basis tetrad 



In the following we shall assume that in each point of space-time, x, there is defined 
a "tetrad", i.e., an orthogonal basis of unit 4-vectors (r, 1, such that the last three 

4-vectors are space-like, and 

e.X^uT'lW'' = 1, (4) 

where e^A^ti/ is the purely antisymmetric tensor. Orientation of the tetrad is arbitrary and 
is defined by 6 sufficiently smooth scalar functions (3 pure space rotations and 3 space- 
time rotations). The other 10 (out of 16 components of the tetrad vectors) fully describe 
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properties of the space-time itself. This feature is used in the so-called "tetrad formalism" 
of general relativity |21|. It is important not to confuse this tetrad with the basis vectors of 
the coordinate system. 

A special choice of orientation links the tetrad (r, 1, 1', 1") to the electromagnetic field 
tensor, F^^. We shall refer to such basis choice as "field- related" . Our subsequent use of the 
tetrad makes other choices possible in the first approximation, jisl but in higher orders just 
the field- related choice of orientation is compatible with the solution, [3] so we shall use it 
from the beginning. 

With the field-related tetrad the (F, l")-plane coincides with the space-like invariant plane 
of the antisymmetric tensor F^^,, while (1, r) result belonging to its other invariant plane, so 
that the Faraday tensor can be fully expressed as jl^ 

F,, = H {I'j; - I'X) + ^ (^/^^^ - ^-^/^) ' (5) 
where H and E are scalar functions, having the physical meaning of the magnetic and electric 
fields in the reference frame where they are parallel. Thus, the field-related basis orientation 
is a generalization of the magnetic coordinates of the non-relativistic treatments. Fradkinjlfil 
first used the invariant planes of the electromagnetic field tensor for decomposition of motion 
of a charged particle. 

The bivector combinations in the expression for F^j^, namely 

can be viewed as linear operators, which combine projection on an invariant plane with an 
orthogonal turn in it. Then the Faraday tensor is 

¥ = Hh + Ec, (7) 

while its dual is 

= -Eh + He. 

The pure projection operators 0^'^\ O^^^ introduced by Fradkin can be expressed as 

C)(") = -bb, O(^) = cc. 

bmce 

C)(«) + = J, and be = 0, it is easy to check that the inverse of the Faraday 
tensor, D'^'^ = {F^^)~^ exists for HE ^ 0, and looks like 

D = -lb + ic. (8) 
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B. Variational principle and the ordering scheme 



Particle dynamics in phase-space can be described on the basis of a variational 
principle. Of course, it is not unique, since different functionals may have coinciding 
extrema. One of the simplest-looking forms of such functional was used by Beklemishev 



and Tessarotto|l3l| for description of particle dynamics in the seven-dimensional phase-space 
in the framework of the general relativity. It allows to describe the particle trajectory in 
the phase-space as a seven-dimensional extremal curve. In the non-relativistic limit the 
extra, seventh dimension is just time, and the whole representation can be easily reduced to 
conventional terms. 

The phase-space trajectory of a particle with the rest-mass and charge in prescribed 
fields can be found from the variational principle 55* = 0, as the extremal of the functional 

S= j g^dx^ = j {qA^ix"") + u,)dx^, (9) 

where q = qa/riiaC^, and is the four- vector potential of the electromagnetic field. Ideally, 
a variational principle in the phase space should yield the standard relativistic equation of 
motion of the particle, plus the relationship between the velocity and displacement, 0. 
As shown in this is indeed the case for the functional (jHl), but only for variations of 
occurring on the hypersurface Uf^u'^ = 1. In this sense the phase space is only seven- 
dimensional. The above variational principle is independent of the parameterization of world 
lines. Furthermore, the above variational principle, as any Lagrangian variational principle, 
is gauge invariant, i.e., any change of the generating differential 1-form of the type 

g^dx'^ ^ g^dx^ + dF, (10) 

where F is an arbitrary smooth function of coordinates and u^, does not change the 
extremal curve. 

In general terms one can define the gyrokinetic transformation as such transformation 
of the phase-space variables, that one of the new variables, called the gyrophase, becomes 
ignorable, i.e., the functional, when expressed in a suitable gauge, becomes independent of it. 
Then, the canonically conjugate variable to the gyrophase, called the "magnetic moment" 
or the "adiabatic invariant" depending on the context, is an integral of motion. As a 
result, the motion of particles in the new gyrokinetic variables is simplified, since the system 
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is effectively integrated in one degree of freedom, eliminating the highest characteristic 
frequency. Once the transformation and equations of motion are found, it is possible to 
simplify the kinetic description as well, i.e., find the gyrokinetic equation (see Sect. VI). 

Construction of the gyrokinetic transformation is done by way of expansion in powers 
of formal small parameters e and A. They are introduced according to the "gyrokinetic" 
ordering scheme , which in our case, following the notation of Ig!], is 

Q^dx^ = {u, + g(^A^(x) + Aa^(x/£))}dx^ (11) 

so that e accounts for the relative strength and inhomogeneity of the background field, while 
A allows distinction between the large-scale background field A^, and the wave- fields given 
by a^. In the final expressions both parameters should be set to 1. 

For an infinitesimal Larmor radius, e, A ^ 0, the particle trajectory is a very tight spiral, 
which allows to represent the transformation as 

x^ = x'^ + Y.e'\^r^,,{y')^ (12) 

s,p 

and 

Ui, = Ui, + w {II COS + sin 0) , (13) 

where 

{u,, w) = «, w') + J2 i^^^ ^)sp ■ (14) 

s,p 

Here = (x'^, 0, /i. My) are the new gyrokinetic variables; x'^ being the coordinate of the 
spiral center, the gyrocenter; the Larmor corrections r^{y'^) describe the difference between 
the particle position and its gyrocenter, and are thus defined to be purely oscillatory in 
the gyrophase, 0; 1' and 1" are the tetrad vectors defining the "plane of rotation", they 
are assumed to vary on the scale of the background field and depend on x' only; is the 
average part of velocity, w' is its "perpendicular" component, while oscillatory corrections 
to the velocity components describe the possibility that the velocity along the Larmor orbit 
may behave in a complex fashion in higher orders. 

The transformation (fT^ . p4|l should preserve the velocity as a unit vector. This require- 
ment, expressed by Eq.Q, is satisfied if Mj, is orthogonal to 1', 1", i.e., 

Uu = UoTy + U\\l^, (15) 
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while parameters w,M||,Mo satisfy 

ul = l + w^ + uii^, (16) 

which follows from u^u^ = 1 for all (j). Note that Uo is not the 0th component of u. 

There are just 6 independent corrections to be defined in each order of the expansion. It 
is possible to introduce the seventh free function by way of a correction to (p, but it is not 
necessary. Thus, the gyrophase (p is defined as an angle in the velocity-subspace , where we 
introduced a cylindrical coordinate system linked to the Faraday tensor. This definition is 
covariant with respect to transformations of the space-time coordinate system, which may 
change the vector components, but not_the vectors themselves. 

As compared to our previous work,[l4j] the gyrokinetic transformation is being searched 
for in broader class of transformations, (|T2|) - |TH) . than before. Our previous formulation 
left components of the velocity unperturbed, i.e., {ui,,w)gp = 0. This made the whole work 
simpler, but there was no solution in one particular limit in orders A*,^ > 1. Without the 
wave-field (A*^) solution existed for all parameters. 

A more natural way of defining the velocity transformation would be to keep Ui, indepen- 
dent of 0, i.e., discard Eq.lfT^. but allow the tetrad to flex with respect to the background 
field tensor. Unfortunately, this involves too many variables, and is not realized at the 
moment. 

The ordering scheme should ensure that all components of the particle velocity are of the 
same order despite the ordering of displacement, ()12|) . This is achieved by ordering the dx^ 
vector components in the new variables as 

dx- = ^d,^=(dx^), + -— d0. (17) 

Here (dx^)^ means the differential of x^ with respect to all variables of the new set except 
0. 

The ordering is designed to ensure that the expansion into the Taylor series in [x^ — x'^), 
i.e., in the Larmor radius, is possible for components of the background field, A^, and the 
orientation vectors, l''^ and l"^. Note, that we are going to expand the vector components 
using partial derivatives, rather than vectors using covariant derivatives. Thus, we are going 
to expand the scalar products of vector fields by basis vectors. For this to be possible, the 
metric coefficients, and, hence, the coordinate system and the curvature of the space-time 
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should be sufficiently smooth, i.e., vary on scales larger than the Larmor radius. In contrast, 
the wave-field, a^(x/£:), is allowed to vary rapidly, so that its Taylor expansion in {x'^ — x'^) 
is impossible. Instead of the Taylor expansion the Fourier representation can be used for 
a^. Nevertheless, can be still Taylor-expanded in powers of the displacement from the 
"unperturbed" orbit. This displacement is proportional to the amplitude of the wave and is 
required to be suitably smaller than the wavelength. 



III. THE GYROKINETIC THEORY 



The goal of the gyrokinetic transformation is to make (j) the ignorable variable, and deter- 
mine the corresponding integral of motion. We shall first suppose that our general template 
of the transformation, ()12|l - (fT7jl . satisfies the requirement, apply it to the variational prin- 
ciple, and then determine the arbitrary functions of the template in such a way that the 
Lagrangian is indeed independent of 0. The procedure will proceed order by order. Most of 
the related algebra is already published, |l| so that we shall drop common expressions and 
highlight differences in the process of derivation. 

First we substitute expressions (fT^ . (|T7jl into (|TT|l and drop terms of order or smaller, 
a.un,.ng that i= of o.det .e.o^nd expanding co,nponent= of the vecto. potential into 
the Taylor series. Then, following |13| we construct the gauge-modified functional 



5G' = Q,dx^ - d Qg(A^ + A;)K + er^)^ 



(18) 



where A'^ = An{x"^). This gauge makes the expression more symmetric and converts all 
references to partial derivatives of the vector potential into components of the Faraday 
tensor at x': 

F (19) 

As a result, 5G' becomes a long expression with terms proportional to dx''^, (dr^)^ , and d0 : 

5G' = (^^A'^ + u,+Xqa, - qrlF,, - eq (^r^^F,,+^r'(rl^^ ^ dx'^+ 

+ e\u^ + Xqa^ - -qrlF^A (drf)^ + e [qK^ + K + ) 

Up to this point there are no differences in the derivation. 
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We proceed by eliminating from 5G' all terms oscillating in 0. This is done by solving for 
displacements r^^ and the velocity corrections {uy.w)^^ in each order in e and A. 

A. Order 

Initially we drop all terms in Eq.()20j). which are of order e^oi higher, while retaining 
contributions of order A. Conditions for to be ignorable for the gauge-modified functional 
5G" = 6G' — edR , where R is an arbitrary Gauge function look like 

M/x + ^l^f^ - Q^iFf,^ = 0; (21) 

|K + A,o„--<F„„)^| =-. (22) 

Here y denotes the oscillating part of namely y = y — y, where y = {y)^ is the 0- 
independent part of y. By a proper choice of the gauge-function R one can always satisfy 
equation (j22I)- This means that 

+ Xqa^ - qrlF^, = (23) 

yields the only essential requirements. The physical meaning of Eq. ()23|) is the relationship 
between the first-order gyro-radius, ri, and the oscillating velocity, ii, along the Larmor 
orbit. The time-like component describes changes in energy due to displacements of charge 
in the electric field. 

If the above requirements ()2H) and ()22|) are satisfied and is ignorable, the hybrid 
variational principle in our approximation can be expressed as 5S" = with 

S" = J 6G" = 1 1 + {u)^ + Xqa,) dx'^ + \{{u, + ^gS,)^^^ d0| . (24) 

Equation can be formally solved for r'( to yield 

r'i = ^D'^^u, + XD'^'-a^, (25) 

where D'^'^ = (F^,^)~^, provided that F^j, is not degenerate, i.e., E ^ 0. 

Unfortunately, the degenerate case of^zero parallel electric field is the only case that is 

in 

studied by all other gyrokinetic theories. ^ Thus we need to present a solution that would be 
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applicable in all cases. To do this we split equation (|23|) in two projections using operators 
b, c . Taking its products with b^^, c^^ we get 

b^^'Uf, + \qh^^~a^ - qHr\^^^ = 0, (26) 
+ + ^^^^'^ = 0' (27) 

where r^^''^ has only (I', l")-components, and r^''^^ - only (1, r). 

We obviously do not need all components of and to satisfy these two equations. So- 
lution is not unique. Our previous choice was to set tt^^"* to zero, which becomes incompatible 
with equations for E = 0. Another possible choice, which is applicable for all regimes, is to 
set r^^^^ = 0. Note, that this choice is one among many with the same applicability, its main 
advantage being just relative simplicity. Different choices lead to non-equivalent gyrokinetic 
transformations (see discussion in Section IV,) in particular, our current choice will make 
the result different from the previous formulation even in the common applicability range. 

From the b-projection, (|^. we immediately get 

= [w{l"^ cos - l'^ sin 0)] ^ + ^b^^a^. (28) 

(here the superscript of ri has been dropped since the other projection, r^'^'', is zero.) The 
c-projection, (j?fj) . yields 



and it follows that 



{u^ + \qa^) = (u^ + \qa^) I'' = 0, (29) 

Uo = -Xqa^T^, 
u\\ = Xqd^l^. 

This means that in orders A° (without the wave) the oscillating parts of w,u\\,Uo can be 
ignored. If A 7^ 0, restriction causes w to have the oscillating part as well. 

After calculating averages with the help of the above solution, the 0-independent func- 
tional 5"' becomes 

S" = l{ (^^; + Ma, + {u)^) dx'^ + /2d0} , (30) 

where 

(u)^ = u\\lfj_ + UoTf, + {w [l'^ COS (p + l'l, sin 0) ) , (31) 
12 



and fi is the wave-field- modified magnetic moment, accurate to order e^, 



w{l'^ COS (p + l^ sin <j)) + ]-Xqa^ 



(32) 



Note tliat only b-projections of the wave-potential enter this definition, unlike the previous 
formulation, where all components were involved. 

The evaluation of the gyrophase averages involving the wave-field is, of course, the 
difficult part. The reason for this is the necessity to transform the given function of space- 
time coordinates a^{x^) into a function of new variables, while the transformation rule 
(123) , (EHl) is itself dependent on a^. Solution of the equation 

= (^'^ + ^ [w{l"^ cos - l'^ sin 0)] ~ + ^h^^~a^ (33) 

can be only obtained by expanding in powers of A. 



B. Order 

Solution of the problem in the limit A = is simple !13|: since the oscillating parts of 
w,M||,Mo can be ignored, 

A^^. (34) 

can be interpreted as the magnetic moment. By varying the functional (jHn|l with respect to 
0, we find d/i = 0, i.e., fi is the first integral of motion in this order. 
Now the functional pup can be rewritten as 

S" = j { (^^; + + n„r^) da;"^ + /id^} , (35) 
with Uo = a/1 + IqHfi + wy^. 

C. Order A^e^ 

This order of the perturbation theory is less trivial, but, due to vanishing w, u\\,Uoi coin- 
cides with the derivation provided in our previous paper, |l|. An important feature arising 
in this order is the solubility condition, 

{ («||/^ + UoT^) D^'^ill COS + C sin 0) }~ = 0, (36) 
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which is satisfied identically if the (!', l")-plane is the invariant plane of F'^'^. This is the 
only requirement that motivates the "field-related" choice of the basis (r, 1, 1', I"). 
In this order the 0-independent functional takes the form 

5" = y I (^^A; + + UoT, + ^efix,^ dx"^ + /id^l , (37) 

where the previous expression for the adiabatic invariant, (jSH), holds. Thus, the second 
order modifications to the Lagrangian can be described as inhomogeneity contributions to 
the effective electromagnetic potential with 

x.-C^-C5;;,-("'« + rr,-^, (38) 

while the form of the adiabatic invariant remains unchanged. Note that some of the effects 
of the gravitational field are also contained in Xti via derivatives of the basis vectors. 



D. Order eU^ 

We proceed by considering the linear approximation in the wave-field amplitude without 
high-order curvature effects. As noted above, the main problem here is the transformation 
of the highly- local wave- field a^(x), where x is the particle position, to the guiding- center 
coordinates (x'", 0, yU, uy). The Taylor-expansion procedure used above for transformation 
of the equilibrium field fails if the wavelength is sufficiently short. Instead, Eq. ()33|) is solved 
using Fourier analysis and expansion in powers of A. 

Assume that the wave field is given in terms of its Fourier components a^(k), then Eq. ()33|l 
becomes 

a^(x) = ^ y"d^A;a^(k)e^'^^'exp|^£A;,(^^(rcos0-rsin0) + Afo^/^fi^^ |, (39) 

where only the second exponential factor depends on 0. 

Here the ordering of the wave- vector k is such that ek ~ 0(1), as usual jsl, but without 
additional restrictions on the "parallel wavelength" and frequency, which appear in tradi- 
tional approaches. Namely, we allow uo/Vt ~ 0(1) and k\\p ~ 0(1) in place of the usual 
ordering ~ 0{e). 

We can calculate the required averages in ()30p and (jH^ by extracting the main exponential 
contribution and expanding it into the Fourier series in cf). We use the following identity 
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exp [i( simjj] = Jn {() Gxp (imjj) , (40) 

n=— oo 

to get 

a^(x) = ——2 / d^A; a^{k)e'^''' ^ J„ (0 exp (m(0o - 0)) exp (^e/^^^r^^) , (41) 

' n=— 00 

where 

Sr" = ^(rcos0-rsm0) +4^"% (42) 
qH H 

characterizes the residual terms. We also defined ^ as 



e = fc±P = {ew/qH) yj{kj"-f + {kj'-f, (43) 



where k± = y {kyl'")^ + (kj'"^)"^ is the length of k_|_, which is the projection of k on the 
( 1', l")-plane; 0o = aictan {kj'"-' / kj'^) is the angle of with respect to the (F, l")-basis, 
p = ew/qH is the Larmor radius. 

To the order e^X^, Sr" in the expression ()4H1 can be set to 0, so that its zero-order 
solution is just 



ar(^) = 7^ / a,ik)e^'^-' ^" (0 exp (m(0o - 0)) • (44) 

' n=— 00 

Note, that neglecting the last exponent in equation ()4H) . or solving by expansion in 
powers of A, is possible only if the particle displacement from the equilibrium orbit due to 
the wave field is much smaller than the wavelength. 

Now that the oscillating part of the wave-potential is defined, we can find the oscillating 
parts of the velocity parameters from Eqs. ()29|l . |T6|) . In particular, in this order 

^(1) = (45) 
w 

with 

u'" = uwl" + UoT". (46) 
Now the averages, present in expressions (jHUj) . become 

«r = i^ / d'ka,{k)e^'^-'jo{0, (47) 

((/'"^sin0 + rcos0)4°^) = ^ jd'^k a,(k)e''^"'iJi(O(/'"sin0o-/""cos0o), (48) 

(*(/;cos0 + /;:sin0)) = -^ J d^A;a.(k)e^^'''ai(O(/;sin0o-^;:cos0o). (49) 
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so that to the order e^X^ 



S" = j { (^^; + Xqa^ + u'^ + u^,^ dx"' + /id0} , (50) 



where = af\ and 



=T?I&M /a^^"' (51) 



2H 



A = - ^^^V.. ) , (52) 



2qH V 

and we have introduced the new notation 

f,-^ = ^J I [k,a, - Ka,] e^kx'M) = ^ J d'k /^,(k)e^^'''^, (53) 

thus expressing the adiabatic invariant fi and uj_^ in terms of the gauge-invariant Faraday 
tensor of the wave, 

ftMu = dai,/dx^ - da^/dx". 

Note that the correction in Eq. ()52|) is proportional to the averaged parallel component 
of the magnetic field of the wave, b'^" f^u- This effect is obviously necessary, but is usually 
neglected due to the standard gyrokinetic ordering of the wave-field, so that the first order 
contribution to the magnetic moment vanishes. 

Note also that in this order we can use the relationship 

ul-w^-ul = l (54) 

instead of (fTBj). 
E. Summary 

Finally, we summarize the above results by presenting the transformed variational prin- 
ciple valid through the first order in A and second order in e, i.e., with terms of the order 
eX and e'^ retained: 6S = yields the particle phase-space trajectory with 

S = j (^qA'^ + u\\l^ + (1 + 2qH*fi + My^) + qa^ + + ^f^X^^ d x'^" + /id0, (55) 

where {x'^,u\\,jl,(f)) or {x''^,u\\,w,(j)) are the new gyrokinetic variables with jj, = w'^/2qH*; 
q = qa/niaC^ is the signed charge-to-mass ratio; 

H* = H + W-f^,, (56) 
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is the effective magnetic field strength, with /^j, being the averaged Faraday tensor of the 
wave, given by Eq. ()53|) . Expressions for the averaged potential of the wave, a^, the wave- 
induced drift, u±fj,, and the second-order correction, x^, are given by equations (jlZ|),(jSH), 
and pHjl respectively. 

The orthogonal basis (r, 1, 1', I") is chosen in such a way that (1', 1") coincides with the 
space-like invariant plane of the antisymmetric tensor of the electromagnetic field F^^ = 
dA'^/dx'^ — dA'^/dx"^, with H being the corresponding eigenvalue. 



F. Equations of motion 

Equations of motion, or the relationships between differentials tangent to the particle 
orbit, can be obtained as Euler equations of the transformed variational principle . Let 
us find the first variation of S assuming (x'", 0, fi, u^\) to be independent: 

SS^ f (-qF^, + d,„)Sx"'dx''' + 

\uo ofi ofi 2 J 

(UwT \ 
+ + s^A [5u\\dx'^ - du\\5x'^] + (5/id0 - d/i(50. 

Here we dropped bars over u\\,uq since all oscillating corrections are left behind, and intro- 
duced a new tensor notation 



(57) 



where the operator Rot,y is defined by Yiotyf^ = df^/dx'^ — dfu/dx'^. Furthermore, the 
partial derivative of the wave-induced drift is denoted as 

_ du^^ _ 1 , - 
'^^ = ^- 2H^^ ^'^^ ■ 

Extrema of the functional are achieved for world lines where (55* = for all variations of 
independent variables. This yields the Euler equations as 

d/i = 0, (58) 
^^-^^^*^)+g?^ + ^dx'-+d0 = 0, (59) 



Mo dfi dp, 

(/m + ^ + dx'^ = 0, (60) 
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{-qF^, + 4m) dx'^ - j^/^ + ^ + s^,j dnii = 0. (61) 

The first equation here confirms that /i is an integral of motion, the second equation deter- 
mines the rate of rotation along the Larmor orbit and is not needed for our purposes here. 
By multiplying the 4-vector equation (jF)T|) by dx'^ and taking the sum in v one can recover 
the third scalar equation ()60|) multiplied by dwy, i.e. out of five equations for four unknown 
functions dx'^ /(\.u\\ only four are independent, as it should be for solvability. This allows to 
determine the non-trivial solutions (at least formally) and use them as coefficients of the 
gyrokinetic equation. 

It is convenient to rewrite the 4-vector equation (jfilj) through its projections on the 
orthogonal vectors of the basis (r, 1, 1', 1"). Taking scalar products with (1', 1") we get 

{qHll + Z'^c/,^) dx'^'+s^.rduw = 0, (62) 
{-qHl'^ + rduf,) dx'^'+s^yrduw = 0. (63) 

These two equations determine the drift velocity. Taking the scalar product with I, we 
recover the following equation governing acceleration parallel to the magnetic field 

dun -qET^.dx" + rdy^dx''' = 0. (64) 

Then, the last equation looks like 

uiidun+UoqElf.dx''' - UoT^di^^dx'^ = 0, (65) 

which is the energy conservation law. However, as shown above, instead one can use another, 
equivalent but much more compact equation (jHOI). It describes the relationship between u\\ 
and the differentials of the guiding-center position 

{uolf, + M||r^ + UoS±f,) dx'^ = 0. (66) 

Equations can be rewritten to exclude dwy using Eq. (jM|) : 

{qHl';^ + rd,^ + s^J" [qEr^ - l-'d,^]) dx'^ = 0, (67) 
{-qHl'^ + rd,^ + s^J'^ [qEr^ - rd.^l) dx'^ = 0. (68) 



18 



It is possible to formally solve equations in terms of the basis vectors. To do 

this we introduce new notations, 



and 



dpq — p di/fiQ^ 1 



Sp = s±^p 



(69) 
(70) 

(71) 



for arbitrary four- vectors p, q. Note that f ^m'^ = and these vectors can be used in place of 



/ r 



Searching for the solution in the form dx'^ =u^ + al'^ + hl"^ + cf^, and solving equations 
(|U m) . (jU7j) .(|UH j) with respect to a, 6, c up to the first order in 1/qH, we get 



qH 



UoV^ 



ul-u\ 



isi"di'u — si'di"u] 



ds', 



(72) 



where ds' is an arbitrary scalar function. Similarly, from Eq. (jf)4|l we find 



du\\ 



{qE - dir) Uods' 
1 



qH 



dw'ditu — dii'di'ti 



UoU\\ 



u„ 



{qE - dir) {si"di>u - si>di>iu) 



ds'. 



(73) 

(74) 



The zero-order terms describe the particle trajectory in quasi-uniform fields. In the non- 
relativistic limit the dir term in Eq. ()73|1 is responsible for the parallel diamagnetic force 
~ fj^VB, while contains the parallel velocity as well as the E x B -drift. The first order 
corrections (in square brackets) describe drifts due to spatial and temporal inhomogeneities. 

Note that the equations of motion in the present formulation result different and signifi- 
cantly more complicated than in the previous formulation. |lj Thus, it should be used only 
if the E = limit is expected to occur in conjunction with the short wavelength wave fields, 
in which case the previous formulation fails. 

In Section V we will need to relate the gyrokinetic displacements to the particle displace- 
ment, i.e., express ds' in terms of the particle displacement ds. Below we present a proof 
that the two displacements coincide in the zero-order approximation. 

According to the definitions. 



ds^ = g^udx^dx" ^ Qi^ui^') 



dx'^ + -^d</. 



dx'^ + ^d0 
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This can be expanded to yield 

d.^ - ,;.dx'-dx- + 2,;.dx-^d0 + ^^d0^. 

The differentials in the right-hand side are supposed to be taken along the particle trajectory, 
i.e., should satisfy the equations of motion, (j3^ . (f72|) . Using them in zero order, 

<d0 + qHr^dx" = 0, (75) 
dx'^" - u'^'ds' = 0, 



we obtain 



q'H'ds'\ 



Furthermore, available expressions for rf can be used to yield 

Or I dri^ 

d(p d(f) q^H"^ ' 

so that, finally, 

ds = ds'. (76) 



IV. EQUIVALENCE OF GYROKINETIC THEORIES 

The gyrokinetic transformation is about finding and using the adiabatic invariant, which 
corresponds to the fast Larmor rotation, and in the gyrokinetic limit becomes the first 
integral of motion. As such it splits the phase space into a sequence of hypersurfaces. The 
invariant hypersurfaces should be the same for all gyrokinetic theories, unless they invoke 
an additional first integral (adiabatic invariant) of some sort. Thus, we conclude, that in 
the same phase space the adiabatic invariants of different gyrokinetic theories should be in 
functional dependence. There is a catch, however. Indeed, the phase space of a canonical 
theory is not gauge-invariant, since the canonical momenta depend on the gauge of the 
electromagnetic field. Thus, the adiabatic invariants with different gauges can (and even 
should) be different. For a non-canonical Lagrangian theory the phase space is chosen in an 
observable way, but this does not preclude similar ambiguity. Indeed, the gauge invariance 
of the Lagrangian itself allows for different expressions to appear. 

There is another sort of ambiguity in terms related to division of the field into the 
"background" and the "wave" parts. The use of the two expansion parameters, e and A, 
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instead of a single one is due to convenience, rather than necessity. It is especially evident 
in our covariant treatment, where all artificial restrictions on the wave- field introduced in 
conventional theories are removed. Indeed, both parameters are inversely proportional to the 
background magnetic field, and the two expansions can be efi^ectively replaced by a single one 
by assuming that all nonuniformity of fields is due to effective "waves", i.e., by representing 
the background field via its Fourier expansion rather than its Taylor expansion. Obviously, 
if the variation of the background field is small on the Larmor scale, i.e., £ <^ 1, then the 
deviation from the unperturbed orbit will also be small, i.e., A ^ 1, and the applicability 
of the wave-field expansion for the variable part of the background field is ensured. Unlike 
the Taylor expansion, Fourier expansion is non-local, and as such has potentially better 
convergence for finite values of the deviation. This is very important for the gyrokinetic 
theory, which is used primarily for description of effects due to the finite Larmor radius of 
particles. Unfortunately, non-local theories are far less convenient from the mathematical 
viewpoint. 

The two ambiguities described above are more-or-less controlled by the user, so that they 
can be counted among advantages rather than difficulties of the theory. If the adiabatic 
invariant is the only integral of motion, and the gyrophase is defined as its canonically 
conjugate angle variable, the freedom of variation becomes confined to different transforma- 
tions within the subspace of the rest of the phase variables. This can be defined as a class 
of equivalence of different gyrokinetic theories: 

//, under the same gauge and the same division of the wave-field, the adiabatic invariants 
of two theories are in functional dependence, then the theories are "equivalent" , 

i.e., the definitions of the gyrophase coincide, and all the difference is just in possible 
transformations of the "drift" variables. As stressed above, theories should be equivalent if 
there are no more integrals of motion that could be expressed via local parameters. (Global 
integrals of motion are here irrelevant, if the derivation procedure is local as usual.) 

The above argument holds for any particular power of expansion in e. However, the defini- 
tion of the adiabatic invariant of order n may depend on the existence of other approximate 
integrals of motion of order n, since any function of all integrals of motion is the integral 
of motion and there is no formal means to distinguish between them. In particular, in the 
first order the electromagnetic field is uniform, and in theories with slow drifts even the 
electric field is zero, so that there are at least two more integrals of motion, such as the 
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parallel velocity and the perpendicular gyrocenter position. As a result, in the first order, 
the definition of the "adiabatic invariant" is very slack, and besides the gauge, allows for 
two more arbitrary functions to enter the definition. This fact is refiected in our paper, jl^ 
where it has been shown that the gyrokinetic transformation in the first order is defined 
with six free parameters. 

For higher-order theories this is dangerous, because by choosing one of the definitions 
at an early stage, one can be left without any solution in the next order, where the other 
approximate integrals of motion may disappear. Thus, the solubility condition for the second 
order approximation removed all freedom from the definition of the adiabatic invariant in 

n n 

our subsequent papers, m, UJ (though left some in other respects.) Another possibility 
is the survival of additional integral (s) of motion in the next order due to some ordering 

n n 

or symmetry assumptions. For example, all previous (except [l|, [IJ]) gyrokinetic theories 
assume that the parallel electric field and the background field inhomogeneity are zero, and 
the parallel wavelength of the perturbation is infinite, while describing the effect of the wave. 
This means, that at least the parallel momentum of the particle is retained as an approximate 
integral of motion in this order, and there is a possible ambiguity in the definition of the 
"adiabatic invariant". (Of course, nobody is so dumb as to insert the parallel momentum in 
the definition by hand. Still, while using multi-step approximations, it is necessary to keep 
the possible ambiguity in mind.) At the moment we do not know of any danger related to 
this possibility, and it seems that, though the adiabatic invariant can be potentially defined 
in different ways, it is impossible to say which one is "right" . Thus, the gyrokinetic theories 
are potentially non-equivalent. In particular, different choices of free parameters in Ref . jl^j] 
lead to non-equivalent theories, when the corresponding definitions of the magnetic moment 
are different. 

So, why is this "non-equivalence" important? The reason lies in the internal machinery of 
the derivation. If the gyrokinetic transformations are non-equivalent, the definitions of the 
gyrophase are different, and, hence, so are the definitions of the gyrocenter/guiding center. 
As a consequence, the "drift trajectories" in different representations will differ by a distance 
of the order of the Larmor radius (xe""^), depending on the gyrophase. Drift trajectories 
are routinely used for confinement analysis. If they form a drift surface, the conclusion is 
that the confinement is achieved. Now we see that this logic works in only one direction. 
The reverse is not true: even if the drift trajectories form a surface in one representation, it 
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will not necessarily be the same one you are using! 

// the gyrokinetic theories are not equivalent, they have close but essentially different drift 
trajectories. 



V. CYCLOTRON DAMPING 

A fundamental issue of for the kinetic description of relativistic magnetized plasmas is 
the consistent treatment of cyclotron losses. In this section we describe the radiation losses 
of a charged particle accelerated by an external electromagnetic field as a single-particle 
classical process. However, one should keep in mind that in cold, dense plasmas there 
may be interference of cyclotron radiation from different particles, thus modifying both the 
outcoming radiation and its net effect on the particle. In this limit the "cyclotron losses" 
are better described by means of collision integrals and plasma waves. 

Radiation damping of the motion of a particle in an electromagnetic field is given by the 



force 
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3c V ds^ ds^ 

Normally it should be a small correction to the electromagnetic force, but can also be large 
for ultra-relativistic particles. 

This expression can be translated to the covariant form by substitution of covariant 
derivatives in place of d^u^/ds^ . Using the equations of motion, we get 

^ 2 7-^ T^kl , D-^ife / k 

-— = q^FikP'^'ui + q—^u'u\ 
ds^ ox'- 

Substituting it into the expression for the force, we can deduce the four-acceleration caused 
by radiation as 



Dm, A 2 



qF^^F^'ui + qu,{F,iu'f + 



(77) 



ds 

where = qi/rnaC^ is the electromagnetic "radius" of the particle, and q = qa/rriaC^ is the 
normalized charge. 

From the normalization it is obvious, that the order of the correction caused by the 
damping is e'^Va/pL, which is, formally, e~^. However, the ratio ra/pi is very small in 
realistic situations, so that we shall assume it sufficient to calculate the trajectory {ui{s)) 
just in the first order in the Larmor-radius expansion. In this order the electromagnetic 
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field, as well as the coordinate system are constant, i.e., the last term in the expression (f77|) 
can be neglected. Thus 

{Du,)^ = '^r^q' - ds. (78) 

This expression has a general form of acceleration due to friction versus some medium 
with characteristic velocity uq = F^u: 

a. = ttr [uq— u (uuo)] . (79) 

It is consistent with ()78|) for a.^ = 2raq^/3, and a= Du/ds. It includes deceleration 
due to both, the cyclotron/synchrotron radiation, which is proportional to H^, and the 
bremsstrahlung, which is proportional to E"^. 

Note also, that the acceleration due to the mean field is proportional to Fu, with F being 
antisymmetric, while the radiative deceleration is proportional to Tu , with a symmetric T, 
but such that uTu= 0. 

Next step is to rewrite the friction in terms of the gyrokinetic variables. Using the 
representations of the world-velocity, 

u„ = u^ + w {ll COS + I'l sin 0) = + u±„, (80) 

and of the Faraday tensor, (|7j), we find 

a = ar [iJ^bbu_|_ + i^^ccuy — u (iJ^u_|_bbu_|_ + i?^U||Ccu||)] . 

Using definitions of b and c, it can be further transformed to 

a = a, [-HW + E\ - u {w^H^ + (1 + w^)E^)] , 

or 

a = -ar {H^ + E^) [{1 + w^)u^_ + w\] . 

Since = —vo^ and u| = 1 + ti?^ it is obvious that au= as desired . 

The cyclotron friction is a small effect, so it can be evaluated in the first-order gyrokinetic 
approximation, i.e., with all fields being uniform. Also, we can assume components of the 
metric tensor to be constant in this approximation. In this way we neglect accelerations due 
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to gravity and the curvature of field lines as compared to the acceleration on the Larmor 
orbit. Then, the equation 

Dur.=ads 



can be rewritten as 



w 



{ll COS0 + /;;sin0) + — 



Ur. 



dWr + 



+ 



1 ^11 
V — ^/^ 

Uo 



&a\\r + {—l'^ sin + I'l cos 0) wdcpr 



-ar {H^ + E"^) [wil + w"^) (/^ COS + /" sin 0) {u\\lf, + UoT^)] ds. 



Taking scalar products with the basis vectors we find 

d0r- = 0, 

dwr = -^Taq^ {H'^ + E^) w{l + w^)ds, 



(81) 



du\\ 



- iraq^ {H^ + E^)w\\ds. 



Here the differentials ds can be identified with ds', which was used in Section III for descrip- 
tion of the gyrokinetic trajectory, due to identity (fT6|l . As a result, it is easy to "tweak" the 
equations of motion to include the radiation damping as a small correction: 



dx'/^ = dx'^ + 
dM|| =du\\g+du\\r 
d/i = + wdwr/qH 
d(j) =d(j)g+d(j)r 



(82) 



where dx^dun 



^j,\\gU.(pg describe the gyrokinetic trajectory, while du\\r,dwr,d(f)r are corrections 
due to the cyclotron damping. All these differentials along the particle orbit are already 
found in Section III and above, and are ready to use. 



VI. KINETIC EQUATION 

If equations of motion contain dissipative terms, i.e., the particle energy and momentum 
are lost along the trajectory, the form of the kinetic equation itself may change. It is 
necessary to check the Liouville's theorem, and introduce corrections if the phase volume is 
not conserved. 
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First, let us define the phase space as a direct product of the space-time by the four- 
velocity space. Note, that scalar products are defined in each subspacc, and between vectors 
of different subspaces, but not in the phase-space as a whole - you cannot add velocities to 
distances. However, all trajectories or world-lines of particles, can be parameterized by the 
proper time s, whose definition involves only the space-time subspace. In effect, zero length 
is ascribed to displacements in velocity directions. 

Coordinates in the phase-space are introduced in accordance with existing scalar prod- 
ucts. Namely, initially, for both subspaces we use the same systems of basis vectors, i.e., for 
an 8-dimensional space we use two copies of four-dimensional basis vectors of space-time. 

Consider a fiux-tube of trajectories. At a given point (x, u) it could be characterized 
by the direction of its central line,(dx,Du)/ds, its cross-section, SV, and the flux density, 
/(x, u), which is identified as the distribution function of particles. However, the definitions 
of the flux density and the cross-section area are subject to the definition of the scalar 
product, which is not fully defined. Ascribing zero length to velocity displacements, as 
in the standard parameterization, means that the direction of the central line is in effect 
(dx, Du)/ds-^ (u,0). Then the phase-space volume 6V can be defined as a product of the 
volume element in the space-time, SV^, which is orthogonal to u, and the element of the total 
velocity subspace, 6Vu, so that SV — SVxSVu. The volume element in the velocity space, 
6Vy,, is also orthogonal to u, since the motion occurs on the hypersurface = 1, which is 
orthogonal to its radius, u. 

The flux of trajectories is continuous, i.e., the flux lines do not end, if the particles do 
not disappear. This means that f{'x,u)6V = const, i.e., 

d(/(x,u)5y) = o, 

which is a general form of the kinetic equation. If, additionally, the Liouville's theorem 
holds, i.e., 

dSV = 0, 

then the kinetic equation obtains its classic form, 

d/(x,u) = 0. 

Note, that the differential d takes place in the 7-dimensional phase-space, while the 
volume element 5V is 6-dimensional, as it describes an element of the hypersurface. Con- 
sequently, / SV is the density of particles in their rest-frame in one particular moment. 
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Defined in this way, the phase-space volume and the distribution function coincide with 
their non-relativistic analogs in the corresponding limit. 

The full volume element in the space-time is 6Va = eaP'ysSxiSx'^SxIdxl, while its sub- 
volume (section ) orthogonal to the unit vector u is 



= ea/B-ysu" Sxf 6x2^x1 + SxqSxJSxI — 6x^6x16x2 — 6xqSx2Sx 



Here 5xf is the /^-component of the i-th box vector, and ea/S'yS is the absolutely antisymmetric 
tensor. Thus, the volume elements are scalars, so that if they are constant in one coordinate 
system, they are constant in all others. In particular, we can choose the local co-moving 
reference frame, in which = and Ua = for all a 7^ (this requires also ^f^o = 0). Then 

= eoi23U° Sx^Sx'^Sx^ = u^^/^5x^5x^5x^ , (83) 

where the box sides are taken along the basis vectors. In words the invariant space subvolume 
is defined as follows: if there is a particle with velocity u, go into its rest frame, and put it 
in a box. This is a unique way of measurement, i.e., a scalar. 

The effective volume element in the velocity-subspace we define as 

5Vu = -e^.^xKu'^Su'^Su^Su^ = -e'"'^^u^SuiJu2x5u3n- (84) 

Its definition is important for calculating moments of the distribution function, such as 
the four-current density. We have defined it as a scalar, i.e., an invariant under Lorentz 
transformations, and it should coincide with corresponding definitions by other authors. 
Indeed, in the co-moving frame it becomes 

= -e''^'uo6u,6u,6u, = ^^^^^^^^ (85) 

rrf-f 

which is listed as an invariant velocity volume in Refs. (l2l . |2l| It can also be expressed via 
the 5-f unction as 

6Vu = 5{l- v/^v^) d^u = 5 (1 - yA^) e'"'^^5uof,5uiJu2x5u3,. (86) 

Indeed, this is also a scalar expression, and in the co-moving frame it obviously coincides 
with Ea. (IHK|l . 

Multiplying the volume elements ()83|1 and (j85|) we get a simple expression for the phase- 
space volume in the co-moving frame 

SV = 5x^6x^6x^6ui6u26u3. (87) 
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Here the box sides Sx"", 5ua are defined along the corresponding basis vectors, i.e., are pure- 
space coordinates and velocities. 

The next question is how this box changes along the particle trajectory. To describe 
it, first note, that each direction in space, n, specifies a rectangle cross-section of the box, 
[xn,Un]. Such cross-scctiou is initially a rectangle, since the volume is defined by a simple 
product of the velocity and space boxes. Consider a small evolutionary change of the box, 
and, hence, of the cross-section. The deformation of the rectangle is small and is due to 
different velocities (in x) and different accelerations (in u) of particles. Since it is small we 
can consider it linear, and as such it is a superposition of deformations proportional to uds 
and to ads = Du. 

The velocity-deformation, x'^ — Xn+Unds, conserves the area, since it acts only in the Xn 
-direction, and u is independent of x. 



Deformations due to perpendicular (to n) components of the velocity cause turns perpendic- 
ular to the plane of the initial rectangle, and thus could affect the area in the second order 
only. 

The acceleration-deformation, u'^ — u^+ands, acts only in the Un-direction, but can 
change the area, since 



Deformations due to perpendicular components of the acceleration have no effect in the first 
order. 

Thus, the change of the area of the cross-section in each direction along the trajectory is 



The change of the volume, 5V, which can be represented as a product of areas of three 
perpendicular cross-sections, is, obviously. 



Note, that 

• due to the choice of the reference frame, the essence of the derivation is tree- 
dimensional, i.e., it is valid in the nonrelativistic case as well; 






d5V = 5Vds ■ divua. 



(88) 
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the Lagrangian phase-space and the volume conservation law, (j89|) . can be generalized 
to multi-particle systems. In this case the volume will be defined as a product of 
volumes for each particle, and the rate of change will be equal to the sum of rates for 
each particle. 

the resulting form can be easily transformed to the general reference frame by replacing 
the three-dimensional divergence by its four- dimensional analog. Adding derivatives 
in the direction of u cannot change the value, since a and u are always perpendicular 
to each other. Thus, 

d5V = SVds ■ Div^a. (89) 

the concept of the four-dimensional divergence is invariant with respect to changes 
from fiat to curved space, and the effect under description is just first order. This 
means, that even if we lost some generality along the way of derivation, expression 
should be valid in the curved space as well; 



• the divergence of acceleration with respect to velocity is very easy to find, since the 
metric tensor is independent of velocity. Indeed, 

• the gravitational (a= 0) and electromagnetic (a" = qF^^^up) forces always satisfy the 
Liouville's theorem, since F°'^gpa = due to antisymmetry of F""^. Also, any force, 
which is independent of velocity, satisfies d6V = automatically. 

The above list of cases, when d6V = 0, i.e., the phase-space volume is conserved, includes 
all most often encountered cases of dynamical systems. Besides, it is still incomplete, but 
also obviously not all-encompassing. In particular, forces due to the cyclotron damping for a 
particle describing the gyrokinetic trajectory cause nonzero contraction of the phase- volume, 

Divua=a, [F'^^Ff^a - 3m„F"^F^^m^] = -a, [H\2 + 3w^) + E\1 + 3w^)] . (91) 

As expected, the phase-space volume of the system decreases due to emission of electro- 
magnetic waves. The rate of the decrease is constant for particles with non-relativistic 
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perpendicular velocities, w ^ 0, but increases as w'^ in the ultra-relativistic limit. The 
kinetic equation for this case looks like 

df + Divua-/ds = 0. (92) 

A. The gyrokinetic equation 

At this moment we are ready to construct the gyrokinetic equation with due account for 
the cyclotron losses. Starting from the kinetic equation for a dissipative system, (jU^ . we 
make a transformation to the gyrokinetic variables, {x''^,u\\fi,(j)), in the phase-space. The 
transformation rules are standard and result in 

df df df df 

■' -dx'^ + TT^dMii + -^d/i + -5-d0 + Div^a-Zds = 0. (93) 



dx'f^ du\\ djl 

At this point we make an assumption that the motion is almost ideal, i.e., the cyclotron 
emission can be described as a small perturbation, and its effects can be calculated in zero 
order in the Larmor radius expansion. As a result, we can define the perturbed trajectory 
and the divergence term via equations (jH^ and (jH^)- A nice common feature of these 
expressions is that they are all independent of the gyrophase 0. This allows to reduce the 
effective number of variables in the following way. 

All coefficients of equation (j93|) are independent of (p, while the solution should be periodic 
in it. If we represent the solution as a Fourier series, different harmonics of the solution 
will be completely independent. Assuming that all 0-dependent harmonics (m > 0) were 
either absent initially, or disappeared due to the phase-mixing, it is possible to consider the 
behavior of the ^-independent part (m = 0) alone. It satisfies, obviously, 

^^-dx''^ + ^d^ii + %dil + Divua-/ds = 0, 



dx'f^ du\\ dp, 

which is the gyrokinetic equation with due account for the cyclotron emission. 

The "magnetic moment" is no longer an exact integral of motion, and is decreasing 
with time. Note that the rate of decrease, and the emission-caused reduction of the phase 
volume produce terms of the same order and are both equally important. In contrast, 
correction to the parallel acceleration seems small and unimportant as compared to its ideal 
counterpart. It is certainly small, but in some cases describes important physical effects. 
For example, if the parallel momentum of the particle is conserved, the cyclotron damping 
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of the perpendicular component of the velocity should be accompanied by the increase of 
its parallel velocity, since the relativistic mass of the particle decreases. 



VII. THE MAXWELL EQUATIONS 

In general, kinetic description of plasmas involves a combination of the kinetic equation 
and the Maxwell equations, which describe the evolution of the collective electromagnetic 
fields. The same is true for the gyrokinetic theory. However, since the gyrokinetic equation is 
written in specific gyrokinetic variables, and thus yields the distribution function in terms of 
these variables, a special procedure of integration in velocity space is necessary to calculate 
the source terms of Maxwell equations. This procedure is presented below. 

The general form of the Maxwell's equations in presence of an arbitrary gravitational 
field is well known |2l|. The second subset 

^ ^ (v^i^^1=-r, (94) 



^—g dx^ c 
has a source term, 

f = cY,(lc^ I ^^fc (x, u) 5K, (95) 

which is the four-current density, expressed via the distribution function of particle species, 
fa, the signed particle charge, qa and the element of the volume of the velocity-subspace, 
5Vu. Our current goal is to express /„ (x, u) in terms of the gyrokinetic distribution function, 
and u^6Vu - in terms of the gyrokinetic variables. 

Starting with 6Vu, we note that the four- velocity can be expressed via Eqs. (fT^ - (fT3jl as 



u,, = w 



[l'^ cos 4> + l'l, sin (f) +u\\lfj, + UoTfj,. (96) 



Taking the box sides along the tetrad vectors, and using Eq.(jH), we find 

, - , wdwd(j)du\\ .. 

SVu = S (v/^ - 1) d^u = (97) 



where Uo = a/1 + w"^ + uy^, and g' = g{x.') is the determinant of the metric tensor in the 
point x', where the tetrad is defined. 

That is all we need in the A = limit, when the differentials in the expression (jHTj) are 
independent. In the presence of a wave, A 7^ 0, the gyrokinetic variables are expressed in 
terms of 10,11^ rather than w,u\\. However, it is more convenient to w,u\^ as independent 
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variablees, but transform the distribution function instead, we can make a transformation 
of variables, wu\i(f) wu\\(f), and apply its Jacobian. In order we retain linear corrections 
only, so that 

wdwdSduw wdivdd)du\\ ( dw du\\ w Un\ , ^ 



Uo Uo \ dw du\\ w Uo 

As a result, the expression for components of the current density can be rewritten as 



Y cqa I [w COS + /"^ sin 0) + + M„r^) (x, u) ^^^^'^^^H , (99) 

^ J Uo 

a 

Further, the distribution function /„ is expressed as the function of the gyrokinetic variables 

fa fa (-^ ^5 A) ^11) 5 

and it is necessary to transform it back to particle coordinates before integrating, as in 
Eq. ()99j) the particle position x, rather than its gyrocenter position x', is kept constant while 
integrating over the particle velocity. This makes it convenient to rewrite Equation ()94|1 as 

- iV^F'^^) = —fV^ = Q^(x), (100) 



dx'^ c 
where the right-hand side is also evaluated at x. Then 

Q^{x) = 471^2 Qa f[w {l'" COS + l'"- sin 0) + u^f + w,r^] ( x - J^eV, J yZ^ ^^^^'^^^H 

a \ i=l / ° 

(101) 

Also, it is necessary to keep in mind that the components of the basis vectors, as well as 
the relationship between the magnetic moment fi and the orbital velocity w, are defined 
through the electromagnetic field tensor at the gyrocenter position, and thus should also be 
transformed back to the particle position via the inverse transformation 

a;/^ = x''' + ^ eVf (102) 

i=l 

before integrating. 

The gyrokinetic transformation has been found by expansion in orders of e and A, and 
this expansion has to be exploited here again. The short-wavelength wave contributions to 
the distribution function can be separated from the slowly changing background by using 
the Fourier components as 

(x"^, A, W||) = /f + fi'^ (k, /i, nil) exp [tk.x"-] d% (103) 
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while Q'^i'x) should be expanded in powers of A and e as 

g'^(x) = ^£U^gj.). (104) 

i,j=Q 

The zero-order current and charge density Q'^qq-^ is easy to find 

Qfoo) = E / H\xn (^l' + rA fi'Hx^,f^,^\)dfcdu\^, (105) 

where the basis components are taken at the current position. Here and below we use the 
following notations: 

111 F)llJ 111^ 

which differ from the magnetic-field invariant due to corrections to the magnetic moment 
(found above) in a suitable approximation. The use of allows to avoid rewriting corre- 
sponding corrections in high-A contributions to the current density. 



A. Order A° 

The background distribution function fa^ is slowly changing on the Larmor-radius / 
gyrotime-scale, and thus can be expanded in the Taylor series around the particle position 

/(°) (.'", i,,«,)= /r /i, «,) - srf ^ - .VJ^ + i-^Vfrrl^ + o(s'), (107) 

while the displacements rf are given by the gyrokinetic transformation as functions of 
(^x'^, ji, u\\) and thus should be expanded as well, 

Here, finally, all functions in the right-hand side can be regarded as functions of the particle 
position x^. The same expansion should be applied to the other factor under the integral, 
namely, 

Gi"^ (v^2^^(/'^cos(/. + r''^sin(/.) +M||P + r^) (109) 

We have 

(a:'^ = (x^ya,«||)-£ri'^+£\^|^ ^-£M^ + ^£VK^^+o(£2), 
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and thus 



er^ — e r 



I 2 U 



a AO) 



^ ^ dx^^ dx^ 



+ G' 



dx'^ 
a2 AO) 
dxf^dx" 



+ 



oie 



The first order displacement of the gyrocenter is given by 

w 

= — (/"'^cos^-r^sin^) + AD^'^a^, 
qli 

and, substituting it into the above expression, we find the first-order current density 

d ///3/2j(0) 



(110) 



(10) 



H 



dx^ 



Similarly, it is necessary to solve the equation |T8|) for ri^ in order to calculate the second- 
order current density. The form of this equation (for A = 0) is such that r2 may have 
components proportional to trigonometric functions of (f) and of 20. However, due to the 
integration in 0, only the first type of terms will produce a non-zero contribution to the 
second-order current density. Thus, it is sufficient to solve for r^i-, which satisfies 



^x'^' 



0, 



or 



1 (JT^ 1 (j 111 

A - -D^A*^.^ = -D^^ {u\\l, + UoT^) ■7^-ij{l""cos(j)- r sin0). 
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q dx'^^ q 

Using this expression we get 



(20) - 



47r^ ttIqC^ J fidjldu 

, , -I 



d 



dx^ 

d //73/2j(0) 



^x'^^ qH 



{ri'^ + ri"') d (^fio)^ 



H 



d ( I 



dx^ 



Ur 



dx"^ \\/H ) dx^ 



(0) 



(112) 



(113) 



(114) 



B. Order A^ 



First we calculate the transform to the particle position variables 



/(I) (x,/i,w||,0) = ^ / f^^ (k,/i,W||) e^'^^exp 



(115) 
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Retaining in the exponent terms of order only [while ke ~ 0(1)], we get 
1 



fi'^ (x) 



47]-2 



w 



-iek. [ -;-(/"^cos(/) - r^sin(/)) + AD'^'^S^ 
' qn 



d^fc. (116) 



The last term in the exponent can be neglected in this order, and the rest can be expanded 
into the Fourier series using Eq. pUl) as 

+00 



/(^) (x) = — / /(^) (k) e"^- Yl (0 e-(^-<^o)d^fc, 

n=— oo 



(117) 



where ^ = k±p, p = ew/qH = e^j2fi/ qH , k± = y {kj"^)'^ + {kJ'"^Y, as before, tan 
{k^l"n I {k/n ■ Then, 



(01) 





1 d^fc/(i) (k) e'^^ 1 


1 Uo . 





l/'^ + MoT^) Jo (0 + 2iefih^''k,^^ 

(118) 

Here we have again used W instead of H to take into account the /i-related corrections in 
future orders. 
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